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Abstract 

We consider the impact of 0(a 2 ) hadronic corrections to the energy spectrum of 
the decay electron in muon decay. We find that the correction can be described, 
within good approximation, by a linear function in the electron energy. Explicit 
expressions for the form factors needed in an approach based on dispersion integrals 
are given. 
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1. Testing the electroweak sector of the Standard Model requires to fix the coupling con- 
stants of the Lagrangian by relating a number of high-precision experiments to theoretical 
predictions. Ideally, experimental data and theoretical predictions should be known with 
comparable precision. Apart from the electromagnetic coupling constant a and the Z- 
boson mass mz, the input of choice is the muon life time r M , with the present best value 
r M = (2.19703 ± 0.00004) fis 0. Experiments are planned at the Paul Scherrer Institute 
|| and the Rutherford Appleton Laboratory || which would reduce the experimental 
error on the muon life time by more than one order of magnitude. 

These experimental data have reached such a precision that quantum corrections can 
be observed. To match this accuracy from the theory side, two-loop radiative correc- 
tions to the muon decay in the full electroweak Standard Model are needed. This is a 
formidable, but not impossible, task. A step in this direction is the calculation of the 
purely electromagnetic corrections to order a 2 in the Fermi theory, which has been per- 
formed by van Ritbergen and Stuart 0, |J (see also f7j). Besides that, also the 0(Nfa 2 ) 
corrections in the Standard Model have been considered in Refs. 0]. 

It is well known that the decoupling theorem is not applicable to the life time calcu- 
lation, which means that it is mandatory to also include the contribution of the heavy 
degrees of freedom. In contrast, the decoupling theorem is fully applicable to the nor- 
malized electron energy spectrum in muon decay, which has also been measured with an 
amazing accuracy at the per-mille level. Again the experimental error is expected to be 
reduced by more than one order of magnitude in a future experiment at TRIUMF 
The spectrum calculation is different from the one for the life time since the Kinoshita- 



Lee-Nauenberg theorem |IJ is not in effect. Consequently, powers of the large logarithm 
hi{m^/m e ) do not cancel in the calculation of the electromagnetic corrections to the spec- 
trum. This becomes obvious when fitting the spectrum corrected to order O(a) to the 
Michel spectrum: the resulting effective Michel parameter differs by about 6 % from its 



lowest-order value [IT], a correction which is more than 10 times larger than the corre- 
sponding correction to the muon life time. At order 0(a 2 ) the radiative corrections may 
be expected to be of the order of several per-mille, i.e. they could possibly be visible in 
present data already, not to speak of future high-precision experiments. 

Given this perspective we present in this note the calculation of the hadronic contribu- 
tion to the energy spectrum of the final-state electron in muon decay. This contribution 
is not expected to be logarithmically enhanced, but nonetheless is required for an even- 
tual complete second-order calculation. The details of the calculation are given in the 
following sections. 

2. We consider the decay of a muon in its rest system, 

// _ (p) -»• e~(p') + v li (q 1 ) + i> e (?2) , (1) 

and define momenta as shown in ([!]). The momentum transferred from the charged 
particles to the neutrino pair, q = q\ + q 2 , is then given by 

q = p-p' . (2) 
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It is convenient to introduce the dimensionless variable 



2E e 

x = 3 

m 

to denote the ratio of the energy of the decay electron E e with respect to the muon mass 
m. Neglecting the electron mass, p' 2 = 0, we see that the kinematically allowed range is 

< x < 1 , (4) 

and one has 

q 2 = (1 — x)m 2 . (5) 

The matrix element Ai for ([!]) in the Fermi theory can be calculated most conveniently 
after a Fierz rearrangement factorizing the amplitude into a current J M which describes 
the [ie transition and a current for the v^Vf, interaction. After squaring and summing 
(averaging) over spins, one can write |.M| 2 as a product of two corresponding tensors. 
The one pertaining to the neutrino interaction can be integrated over the unobserved 
momenta of the neutrinos independently, leading to 

N^v = qtfu - 9^q 2 • (6) 

This tensor will be contracted with 

Cpv = J*J U , with J M = u e (p')A M (g)w M (p) , (7) 

where A M (g) is the effective vertex of the four-fermion interaction. At the lowest order, 
A At (g) is identified with 

K = ^7, (1 - 75) , (8) 

where Gp = (1.16637 ± 0.00001) x 10 -5 GeV~ 2 is the Fermi coupling constant. 

Hadronic contributions to the radiative corrections to the current J M at order 0(a 2 ) are 
described by the Feynman diagrams shown in Fig. [I]. The hadronic vacuum polarization 



is inserted in a one-loop vertex correction. The vacuum polarization can be related to the 
measured cross section for e + e~ — > hadrons with the help of a dispersion relation 

U^(k 2 ) = ^ T ^R(s) k \ (10) 
Sir Js thT s k z — s + lO 

where 

<j(s; e + e~ — ► hadrons) 
R ( s ) = — 7 + -\ 11 
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and the integration starts at the two-pion threshold, Sthr = 4m^. Therefore, the calcula- 
tion corresponds to a one-loop vertex correction with a photon of mass y/s, i.e. using a 
propagator 

" V - • (12) 



k 2 - s + iO V " k 2 , 
The result can be written in the form 

/y POO Aq _ 

A M (g) = -/ -i2( S )A M ( S ; g ,m 2 ) (13) 
where the vector function A M can be decomposed into Lorentz-covariants as 



1 I 2\ [-1 | 7 ! 2 2\1 , ~ / 2 2\ , ?C ~ / 2 2\ 

A M (s; q, m ) = ~i^u L 1 + f{s;q ,m ) H -u R g + (s; q ,m ) -\ u R g_{s; q ,m ) , 

(14) 

with = (1 ±7s)/2. The calculation is straightforward and corresponds to that of the 
one- loop vertex correction in QED, with the difference that (i) the exchanged "photon" 
is massive with mass ^/s, (ii) the coupling is purely left-handed, and (Hi) the two fermion 
lines have different masses. In fact, except for small q 2 the electron mass can safely be 
neglected. Explicit results for the form factors / and g± are given in the appendix. After 
contraction with the neutrino tensor N^, Eq. only the form factors f(s; q 2 ,m 2 ) and 
g+(s; q 2 ,m 2 ) will remain in the final result. 

Performing the dispersion integral one obtains the form factors 
f(q 2 ,m 2 ) = ^- [" ^R(s)f(s;q*,m 2 ), 

07T Js thv S (15) 

g + (q,m) = — —R{s)g + {s;q,m). 

07T J Sthr ^ 

Since we are interested in the 0(a 2 ) correction, it is sufficient to keep only terms of first 
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order in / and g + in the decay spectrum which then can be written in the form 



^dT 

r dx 



2x 2 [(3 - 2x) (1 + 2f(x)) + xg + (x)} 

i dr 



r dx 



[l + r(x)) 



(16) 



with 



and 



r r 



G 2 F m 5 

1927T 3 



Born 

dr 

dx 



Born 



r (x) = 2f(x) + 



x 



2x 2 (3 - 2x)r c 



(17) 



where f{x) 



3-2z* 

/ ((1 — x)m 2 , m?) and g+(x) = g+ ((1 — x)m 2 , m 2 ). 



3. For our purpose, the function R(s) describing the hadronic cross section of e + e~ anni- 
hilation can be modeled by a combination of experimental data and analytical results from 
perturbative QCD. Since we are going to calculate a small correction, it is not necessary 
to invoke the most sophisticated treatment as needed, for example, when calculating the 
hadronic contribution to the fine structure constant a(m^). At low s < 2.5 GeV 2 we use 
experimental data from ALEPH parametrized in jjHJ] or provided directly by ALEPH 
from a measurement of the isovector r spectral function. These data are complemented 
by the resonance contributions from the isospin-0 light mesons u and 0. Above s = 2.5 
GeV 2 we use the QCD prediction for R(s) due to light quarks at order 0(a s ). 

Since the data in the cc-channel published by various groups are in a large part of 
the energy range inconsistent, we apply in this case the QCD-based approach of analytic 
continuation by duality [14| . The data region can be chosen to extend only over the sub- 
threshold resonances, i.e. one can calculate the contribution coming from the cc-channel by 
a combination of data describing the J/^f(lS) and J/^f(2S) resonances and the prediction 
of perturbative QCD. We checked that the results obtained this way are consistent with 



those of the standard approach using the new BES data [R 
The correction to the total decay rate, 

Jo dx 

was calculated before in [Q. Our result, 



(19) 



Ar had ~ -0.0421 - r , 



a 



71 



(20) 



agrees perfectly with the corresponding number —0.042 given in ||. The resulting correc- 
tions to the spectrum (16) are shown in Fig. ^|. At small x, the corrections are positive, 
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Figure 2: Results for the form factors f , g + and r defined in (\Th]), flfify. 

but the correction to the total decay width is dominated by the negative values at x^O.18. 
The dependence of the form factors on a; is to a very good approximation linear: 



For g + , the coefficient of the term linear in x is very small and is therefore omitted. Note 
that this behaviour cannot be described by a simple redefinition of the Michel parameter. 
Since Gp is a free parameter in the Fermi theory, the correction to the total decay width 
is not observable; it can be absorbed by a suitable redefinition of the Fermi constant. 
However, the modification of the spectrum is, in principle, observable. 

The correction to the total decay rate can be split up into the various contributions to 
the hadronic vacuum polarization, as shown in Table [1|. The form factor / of the 7 M term 
contributes —0.0387, whereas the correction due to g + , —0.0033, is smaller by one order of 
magnitude. The total correction (the contributions due to / and g + ) is saturated to 81 % 
(80.8% and 89%, respectively) by the contributions from small s below 2.5 GeV 2 . Only 
5 % of AFhad is due to charmed states, and the bottom sector is completely negligible. 




(21) 
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Contributions to Al^ad 



1 

1 


< s < 0.2 (orev 


n nn ion 
— 0.00129 


010/ 
3.1 7o 


2 




-0.00223 


5.3 % 


3 




-0.00264 


6.3 % 


4 


0.2 < s < 2.5 GeV 2 


-0.02804 


66.6 % 


5 


s > 2.5 GeV 2 


-0.00564 


13.4 % 


6 


J/*(15) 


-0.00066 


1.6 % 


7 


J/V(2S) 


-0.00017 


0.4 % 


8 


charm, s > 4m 2 


-0.00138 


3.3 % 


9 


bottom, s > 4m 2 


-0.00003 


0.1 % 




Sum 


-0.04207 


100 % 



Table 1: Contributions to the corrections of the total decay rate 



The numerical results given here take into account the 0(a s ) QCD corrections in R(s). 
Using the leading-order expression for R(s), i.e. without the correction factor (l + a s /7r) in 
the light-quark contribution at large s and the charm-quark contribution, the final result 
would be —0.04149, i.e. changed by 1.4%. We conclude that a more refined treatment 
which would include higher orders of perturbative QCD and mass-dependent corrections 
in the heavy-quark sector is not required for our purpose. 

The evaluation of the dispersion integrals has been performed using standard numerical 
integration routines up to a value s max of several hundred GeV 2 . The contribution above 
this value was obtained with the help of Maple using the asymptotic expansion of the 
form factors for large s (see Eqs. (35)-fl3?D in the Appendix). For intermediate values of 



s, good consistency of both procedures has been verified. 

4. The same set of formulae can be used to calculate the contributions from a / u + 'yT or 
a t + t - loop insertion. We find that the tau loop gives a very small contribution, about 
1.5 % of the one from the muon loop, in agreement with Ref. ||. Therefore we give only 
results for the muon loop, where one has to insert in Eq. (|T3j) 

and s t hr - > 4m 2 . (22) 

V * J ¥ * 

With this input we obtain 

Ar muon ~ -0.0364 (-) 2 T , (23) 




TT 



which perfectly agrees with the exact result given in Ref. |J. The results of a linear fit of 
the form factors for the muon-loop insertion are: 



Uuon(x) ~ (0.0130- 0.0414a;) J T 



2 



7 



9+, 



x 



-0.0090 + 



.ooo5x) 



(24) 



(0.0267 - 0.0898x) ( ^ ) T 



a: 



5. To summarize, we found that the energy spectrum in the decay of an unpolarized muon 
is corrected by a smooth function in x due to hadronic contributions at order 0(a 2 ). At 
both ends of the spectrum no particularly outstanding enhancement or suppression is 
observed. 



The calculations described in this paper constitute only the most straightforward part 
of a full calculation which would be necessary before the expected future high-precision 
data can be confronted with theoretical predictions. This will not only be necessary for a 
meaningful test of the electroweak Standard Model, but also when searching for physics 



beyond the Standard Model [TH . 



Appendix 

The form factors introduced in (|14]) are given explicitly in the following: 



f(s; q 2 ,m 2 ^ 



- — { 2(m 2 -q 2 -s) 
47T v ' 



1 + 



q 2 s 



[m 2 — q 2 ) 2 



Co 



+ 



2q 2 s 



2q 2 



(m 2 — q 2 ) 2 (m 2 — q 2 ) 
2q 2 s{m 2 + q 2 ) 

+7sE(s;g 2 ,m 2 ) , 



+ 1 



B (q 2 ;m 2 ,Q) 



+ 



B (m 2 ;m 2 ,s) + 



m 2 — q 2 



A(s) + 2 



(25) 



g+(s;q ,m ) 



a m 2 ( 2q 2 s 
An q 2 1 m? — q 2 



3q 2 s 



(m 2 — q 2 ) 2 



+ 2 



(f 



m 2 — q 2 



m 2 — q 2 



m 2 — q 2 



6q 2 s 



{m 2 — q 2 ) 2 



B (q 2 ;m 2 : Q) 



m 2 — q 2 



6m 2 s s(4m 2 — q 2 



{m 2 — q 2 ) 2 m 2 {m 2 — q 2 ) 



A(s) + 



+ 2 



m 2 — q 2 in 2 



B (m 2 ]m 2 ,s) 
' A(m 2 



m 2 — q 2 m 2 



(26) 



8 



,2 2 



g-[s;q ,m ) = - T 



a m 2 \ 2q 2 s 



An q' 2 I m 2 — q 2 



3q 2 s 



2s 



(m 2 — q 2 ) 2 m 2 — q 2 



+ 2 



Cn- 



m 2 — q 2 



1 



m 2 — q 2 



m 2 — q 2 



6m 2 q 2 s 2q z s _ 2 „ 2 
H ^ o + 2m 2 - 3q 2 



(m 2 — q 2 ) 2 m 2 — q 2 



6m 2 s s 
+ 



B (q 2 ;m 2 ,0) 



-| h 2 

(m 2 — g 2 ) 2 m 2 — q 2 m 2 



B (m 2 ;m 2 ,s) 



1 



m 2 — q 2 \m 2 m 2 — q 2 



Ms) + 



m 2 m 2 — g 2 



A(m 2 ) , (27) 



where Co = C*o(m 2 , 0, g 2 ; m 2 , s, 0) is the three-point integral (cf. Fig. [I|a), defined in 
Eq. (|34D below. _B and A denote the tadpole and two-point integrals |17|, [18j , respectively 
(see Eqs. ©-© below). 



Self-energy diagrams (cf. Fig. |T]b, c) contribute to the coefficient of 7^ only and are 
given by 



a 







f S E(s;q\m 2 ) = 2( S + 2m 2 ) — B (p 2 ;m 2 } s) 



57T 



dp 2 



2 2 m 



^B (m 2 ;m 2 } s) 



2 1 3" 

-A(s) + — A(m 2 ) + 1 + - 



(2f 



where the last term, |, comes from the self energy on the massless (electron) leg. Using 
recurrence relations |ff9| (see also Appendix A of pOfl), the derivative in (|28| ) can be 
represented as 



_d_ 

dp* 



B (p 2 ;m 2 ,s) 



m 2 (s — Am 2 ) 
+A(m 2 ) - I ] 



-(s-3m 2 )B (m 2 ;m 2 ,s) 
A(s) + m 2 



2m 2 * 



The required tadpole and two-point integrals are JT7|, |1"8 

,2 



Aim 2 



m 



-A - 1 + In 



m 



B (m 2 ;m 2 1 s) = A + 2 - In 



m 



S 1 S S n 1 

In — + — —ft In 



'1 + Ps 



yUj5 R 2m 2 m 2 2m 2 y 1 — ft 



,2 „2 



5 (g 2 ;m 2 ,0) = A + 2 - In -=- + ^ln 

A^dr 9 2 



,2 „2 S 



where 



Am 2 



(29) 



(30) 
(31) 
(32) 

(33) 



and /iDR is the scale parameter of dimensional regularization. In Eqs. (|30"D-(|32j), terms 
containing A = 1/e — In 7r — 7e represent the ultraviolet singularities which cancel in the 
final results (E5D-(E7D. 



Finally, we need the three-point scalar function Cq |17j for positive values of q 2 . In 
this case it can be written in the following form: 

,2\ / ( m 2 _ g 2 ) 2 \ 



C (m 2 ,0,g 2 ;m 2 ,s,0) 



1 

m 2 — q 2 



Li 2 1 



m 



Q 



Li 9 



sq z 



-Li 2 



+- In 
2 



+ - In 

2 



Li 2 



(m 2 — g 2 ) 2 
sm 2 



In 



1 + 



(m 2 — g 2 ) 2 



m 2 + g 2 + (g 2 — m 2 )j3 s 
m 2 + g 2 — (g 2 — m 2 )(5 s 



sq z 



In 



i 2 2 

1 m m 
In — - In — 



2 g 2 



'1 + P, 
I -13s 



(34) 



For convenience, we also give asymptotic expansions of the form factors (|25|)-(|2T|) valid 
for large s, 



4 -f 

a 



3s(q 2 ) 2 

2 



(m 2 -g 2 ) 2 (m 2 + 2g 2 )ln 



' m 2 — q 2 
m 2 



(g 2 ) 2 (3m 2 -2g 2 )ln^- 



rrr 



+q 

4 



4 ^ 2 2 H / 2\2 

m m q [q 

2 H 3 V . 



1 



6s 2 (g 2 ) 3 

2 



(m 2 -g 2 ) 4 (m 2 +g 2 )ln f ^-^-\-{ q 2 f [ 3 5m 4 -3m 2 g 2 + (g 2 ) 2 



In 



+9 



m 8 - ^m 6 g 2 + ^m 4 (g 2 ) 2 + ^m 2 (g 2 ) 3 - ^(g 2 ) 4 



>+0(s- 



(35) 



47T „ 

9-\ 

a 



m 



m 



3s Qs 2 (q 2 Y 



(m 2 - g 2 ) 4 In ( H^l) + ( g 2 )3 (4m 2 _ g2) ln « 



m 6_Z 4 2_8 2( 2)2 _5_ ( 2)3 

2 3 v ; 12 v ; . 



(36) 



— 



4m 2 



,2 „2 N 



3s(g 2 )' 



(m 2 -g 2 ) 3 ln' m * 1 ' 2 



+ (g 2 ) 3 ln — + g 2 



4 ^ 2 2,'/ 2\2 

m --m g +y^(<? ) 



6s 2 (g 2 ) 4 



(m 2 - g 2 ) 4 (6m 2 - 5g 2 ) ln ( — / ) - (g 2 ) 4 (26m 2 - 5g 2 ) ln ■ * 



6m 8 - 26m 6 g 2 + -m 4 (g 2 ) 2 - -m 2 (g 2 ) 3 + -(g 2 ) 4 
2 6 12 



>+C(s" 



(37) 
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These expressions turned out to be useful for the numerical evaluation of the dispersion 
integrals in the large-s region. 

The asymptotic values for x = are as follows: 



a 



' ' l s + 2m 2 )(7s-22m 2 )^- In ( J + (18m 4 -6m 2 s-7s 2 ) In Aj 

x =o rzin"* l p s \1 — B S J m 2 



47T 

—9+ 

a 



47T 

—fi- 
at 



x=0 



+m 2 (33m 2 + 14s) 



1 h{2m 2 -s)f3M ) + (2m 4 -4m 2 s+s 2 )ln^- + m 2 (5m 2 -2s) 



4m 4 1 \1—(3 S I m 



' > s(2m 2 - 5s)(3 s \n I 1 + - (6m 4 + 12m 2 s - 5s 2 ) In 



.=0 12m 4 { v V-Ps) m 2 

+m 2 (9m 2 - 10s)|. 

Note that there are no contributions of the order 0(xlnx). 
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